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Abstract 

One of the challenges in quantum information is the 
demonstration of quantum coherence in the opera- 
tions of experimental devices. While full quantum 
process tomography can do the job, it is both cum- 
bersome and unintuitive. In this presentation, we 
show that a surprisingly detailed and intuitively ac- 
cessible characterization of errors is possible by mea- 
suring the error statistics of only two complementary 
classical operations of a quantum gate. 

1 Introduction 

Due to rapid technological advances, there is now 
an increasing number of experimental realizations of 
multi qubit quantum processes. As more and more 
complex processes become feasible, it is desirable to 
develop efficient strategies for the characterization of 
such gates. Specifically, it may be useful to identify 
the noise characteristics of a gate without having to 
perform the large number of measurements necessary 
for full quantum process tomography. 

A particularly promising approach appears to be 
the analysis of a pair of complementary classical op- 
erations performed by the quantum device 01 El El- 
The errors observed in such a pair of classical oper- 
ations provide a surprisingly detailed "map" of the 
noise in the complete quantum operation. It is there- 
fore possible to use this limited set of measurement 
results to estimate the performance of the gate in op- 
erations that have not been evaluated explicitly. 

In the following, we will briefly explain why the 
coherence of a quantum operation is completely de- 
scribed by the "parallel" performance of a pair of 
complementary operations, and how the observable 
error syndromes correspond to the process matrix ele- 
ments describing the complete quantum process. We 
then apply our analysis to the actual experimental 



data obtained from a recently realized optical quan- 
tum controlled-NOT gate. In particular, we can de- 
rive estimates for the success of entanglement gen- 
eration and of Bell state analysis, even though the 
operations actually measured were completely local. 

2 Error classification 

In general, any noisy quantum process in a d dimen- 
sional Hilbert space can be described using a set of 
d 2 orthogonal operators Aj by treating the operators 
as vectors in an operator space with an inner prod- 
uct defined by the product trace of two operators. 
If the intended operation is given by the operator 
Uq, errors may be classified according to deviations 
from the intended output by multiplying an orthog- 
onal set of output error operators Aj with the ideal 
operation Uq. A specific process E is then described 
by a process matrix with elements so that the 
input-output relation is given by 

Pout = E (Pm) = Xij&iftoPirfto&i- C 1 ) 

id 

In principle, any choice of operator basis Aj is al- 
lowed. For practical purposes, however, Aj should 
be as close to the directly observed errors as pos- 
sible. For operations on quantum bits, the process 
should therefore be expanded in terms of the error 
syndromes given by the Pauli matrices /, X, Y, and 
Z, acting independently on each qubit. Since quan- 
tum information processes are usually defined in the 
computational basis given by the Z eigenstates (the 
Z basis, for short), it is possible to identify X with 
a bit flip error, Z with a phase error, and Y = iXZ 
with a combined bit flip and phase error [5]. How- 
ever, the symmetry of the operations suggests that 
it may be more realistic to interpret each operator 
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as a Bloch vector rotation by tt around the appropri- 
ate axis. It is then clear that a phase error in the Z 
basis will show up as a bit flip in the X basis, ans 
vice versa. Since experiments on individual qubits 
can only measure the bit value of a given basis, it 
is more realistic to regard "phase errors" as observ- 
able bit flips in the complementary basis. It is then 
possible to identify each error operator by a pair of 
observable error patterns, one in the computational 
Z basis and another in the complementary Z basis. 



Table 1: Examples of error indices based on the error 
patterns observed in the complementary Z and X 
basis. 



Error Aj 


Z error 


fx 


X error f x 


index i 


I®X(g>Y 


Oil = 


3 


001 = 1 


3,1 


Z®Y ®Y 


Oil = 


3 


111 = 7 


3,7 


Y %I®Y 


101 = 


5 


101 = 5 


5,5 



A complete set of d 2 = 2 2N orthogonal errors Aj is 
thus obtained by simply combining the 2^ possible 
N qubit errors in Z with the 2 N possible errors in 
X . For convenience, the index i can then be written 
as i = (f z , f x ), where f z / x is the value of the binary 
number obtained by assigning a digit of 1 to the loca- 
tions of bit flip errors. Tabled shows some examples 
of this labeling system for three qubit errors (64 pos- 
sibilities). 



3 Evaluation of experimental 
results 

An experimental quantum gate can now be tested by 
using only two complementary sets of inputs, {| n}} 
and {| A;}}, chosen in such a way that the correct 
outputs of the d = 2 N orthogonal states | n) are the 
cigcnstatcs of Z, {Uo | n) =\ Z n )}, and the correct 
outputs of the d — 2 N orthogonal states | fc) are the 
eigenstates of X, {U n \ k) —\ Xk)}. The experimen- 
tal results can be represented in a pair of error tables 
showing the probabilities for the various outcomes. 
Each outcome can then be identified by its error num- 
ber f z / x , so that the correct output Z n (Xk) is given 
by the conditional probability p(0\Z n ) (p(0\Xk)), and 
the observed error probabilities are given by 

p{f z \Z n ) = (Z n | k fmt0 E(\ n)(n |)A /a , | Z n ) 
p(f x \X k ) = (X k \A 0Jx E(\k)(k\)A 0Jx \X k ). 



Outcomes corresponding to the same kind of error 
can then be averaged to obtain the fidelities Fz and 
Fx and the error probabilities rjz(fz) and T)x{fx) of 
the Z and X operations, respectively. Specifically, 
the sums defining the fidelities and the error proba- 
bilities read 



^ d-l 
Fz= d^ P{0lZn) 
d-l 



d-l 



Vz{fz) = \ E P(f^ 
z n =o 

d-l 

Vx(fx) = -,J2 P(f*\ X k)- 



x k =o 



x k =o 



(3) 



Table [3 shows the experimental data reported for 
an optical quantum controlled-NOT 2 arranged ac- 
cording to the errors f z and f x . Note that the perfor- 
mance of the experimental device is now described by 
only two fidelities and 2d — 2 error probabilities. For 
the quantum controlled-NOT and other two qubit op- 
erations, this means that only 8 characteristic proba- 
bilities are used to evaluate a process fully described 
by a total of 256 process matrix elements. 

Classical intuition already indicates that the prob- 
abilities of the errors f z observed in Z should corre- 
spond to the sums of the diagonal process matrix ele- 
ments with the same value of f z in the first part of the 
index i, and the probabilities observed in X should 
correspond to the sums diagonal process matrix ele- 
ments with the same f x . This relation can indeed be 
confirmed by applying the process matrix definition 
in eq.Q to the averages of the conditional probabil- 
ities defined by eqs. and (0. As expected, the 
results read 



d-l 



d-l 



(2) 



p z = E X(o,/*)(o,/»)i Vz(fz) - E XU*,f»)U*,f*)> 

d-l d-l 

F X = E X(/*.0)(/*,0)) Vx(fx) = E X(f*,M(f*,f*)- 
/*=o /,=0 

(4) 

It is therefore possible to interpret the diagonal el- 
ements X(f z ,f x )(f z ,f x ) as joint probabilities of the er- 
rors f x and f z , where the unconditional probabilities 
of f z are given by the measurement results Fz and 
f)z(fz), and the unconditional probabilities of f x are 
given by the measurement results Fx and rjx(fx)- 
The process matrix elements can then be estimated 
by making some additional assumptions about the 
correlations of errors in X and in Z. 

The relation between the observed error distribu- 
tions and the diagonal elements of the density matrix 
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Table 2: Measurement data for an experimental two qubit gate 2 arranged according to the observed output 
errors. The error distribution is summarized by the averages at the bottom of each table. 



P{fz\Z n ) 


fz 


= 


fz = l 


fz = 2 


fz = 3 


z n = oo 


0.<< 


398 


0.031 


0.061 


0.011 


z n = oi 


0.<< 


385 


0.021 


0.088 


0.006 


Z n = 10 


0.< l 


319 


0.054 


0.031 


0.096 


z n = u 


0.< l 


310 


0.099 


0.027 


0.064 


averages 


F z = 


0.853 


773(1) =0.051 


r) Z (2) = 0.052 


Vz (3) = 0.044 


P(fx\X n ) 


fx 


= 


fx = l 


fx = 2 


/x = 3 


x k = oo 


OA 


354 


0.044 


0.063 


0.039 


x k = oi 


0.1 


370 


0.019 


0.071 


0.040 


x k = 10 


0.1 


371 


0.058 


0.050 


0.021 


x k = u 


0.1 


374 


0.013 


0.099 


0.013 


averages 


F x = 


0.867 


77^(1) = 0.034 


Vx (2) = 0.071 


r) x (3) = 0.028 



can be visualized by arranging the diagonal elements 
X(f,J m )(f,j x ) m a table so that the lines represent 
the errors f z observed in Z and the columns repre- 
sent the errors f x observed in X. Table |3] shows this 
arrangement for a two qubit operation such as the 
quantum controllcd-NOT that was first analyzed by 
this method in |2|. The experimentally observed fi- 
delities and error probabilities are then given by the 
sums of the corresponding column or line. 

Since all diagonal elements of the process matrix 
must be positive, the measurement values impose 
rather strict limitations on various properties of the 
gate. Of particular interest the estimate of the pro- 
cess fidelity X(o.o)(o.o) = Fqp-i since it provides a mea- 
sure of how close the experimental gate is to the in- 
tended ideal process defined by Uq. The relation be- 
tween the process fidelity F qp and the experimentally 
observed complementary fidelities Fz and Fx can be 
determined directly from eq.lQJ, making use of the 
fact that the trace of the process matrix is one. The 
sum of the two complementary fidelities then reads 

(d-i d-i \ 
E E *(/.,/. x/../-) • 
/»=i/*=i / 

(5) 

Since the sum over process matrix elements with both 
/„ 7^ and f x ^ is always positive, the sum of Fz 
and Fx can never be larger than the process fidelity 
plus one. At the same time, eq.Q ensures that no 



classical fidelity is smaller than the process fidelity. 
For this reason, the two measurement results Fz and 
Fx limit the process fidelity of an experimental de- 
vice to the interval given by yQ 

F z + F x - 1 < F gp < Mm{F z , F x }. (6) 

The difference between the upper and the lower limit 
of this interval is equal to Max{F^, Fx} — 1. For de- 
vices with at least one very high fidelity, this estimate 
is therefore nearly as good as a direct measurement 
of the process fidelity. 

In principle, the process fidelity can then be used 
to estimate important gate properties such as the en- 
tanglement capability of the gate 0- However, more 
precise information is available if the complete error 
statistics are used to estimate the fidelities of oper- 
ations other than the Z and X operations observed 
in the experiment [3JE]. For this purpose, it is use- 
ful to construct a "worst case" error model, where all 
errors are either pure Z errors (i = (/ z ,0)) or pure 
X-errors (i = (0, f x ))- The diagonal elements of the 
process matrix are then given by the corresponding 
measurement results, 

X(/„o)(/„o) = Vz(fz), 

X(o,/ x )(o,/ x ) = Vx(fx), (7) 

and the process fidelity has its minimal value of F qp = 
Fz + Fx — f ■ In the case of the quantum controlled- 
NOT gate |2 , the process matrix elements defined by 
this noise model are shown in table 0] 
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Table 3: Illustration of the limits on diagonal process matrix elements denned by the measurement results 
for the complementary Z and X operations for a general 2 qubit operation. 



XCf. ,/»)(/«>/«) 


/x=0 


/x = l 


fx = 2 


/* = 3 


Sum 


fz = 

/. = 1 

/* = 2 
/, = 3 


X(i,o)(i,o) 

X(2,0)(2,0) 
X(3,0)(3,0) 


X(o,i)(o,i) 
X(i,D(i,i) 

X(2,l)(2,l) 
X(3,l)(3,l) 


X(0,2)(0,2) 
X(l,2)(l,2) 
X(2,2)(2,2) 
X(3,2)(3,2) 


X(0,3)(0,3) 
X(l,3)(l,3) 
X(2,3)(2,3) 
X(3,3)(3,3) 


^(1) 
^(2) 
Vz(S) 


Sum 


Fx 




Vx(2) 


^x(3) 


1 



4 Noise models and fidelity es- 
timates 

It is now possible to make predictions for other op- 
erations based on this noise model. In particular, 
any operation resulting in local output states that 
are eigenstates of some combination of X, Y, and Z 
cigenstates has a fidelity given by a well defined sum 
of d = 2 N diagonal elements of the process matrix, 
corresponding to the error operators that stabilize the 
output states. For example, the two qubit operation 
resulting in ZX outputs has a fidelity of 

F zx = F qp + X(1,0)(1,0) + X(0,2)(0,2) + X(l,2)(l,2) (8) 

In the case of the quantum controlled-NOT, the 
ZX operation is simply the identity operation, since 
the ZX eigenstates are also eigenstates of the ideal 
controlled-NOT operation Uq. We can therefore esti- 
mate how well the quantum gate preserves its eigen- 
states. In terms of the error probabilities r] x / z (F x / z ), 
the result reads 

F zx > l-T 1z (2)-r ) z(3)-r ]x (l)-r lx (3) = 0.842. (9) 

Note that this result is significantly larger than the 
result of 0.72 defined by the minimal process fidelity 
Another operation of great interest is the genera- 
tion of entanglement from local XZ inputs. The ideal 
operation Uo converts these input states into the four 
orthogonal Bell states, characterized by their XX, 
YY, and ZZ correlations. The fidelity F^ of this 
operation is also given by a sum of four process matrix 
elements, corresponding to the Bell state stabilizers 
II, XX, YY, and ZZ, 

F El = F 1P + X(3,0)(3,0) + X(0,3)(0,3) + X(3,3)(3,3)- (10) 

The measurement results define a minimal fidelity for 
this operation given by one minus the single qubit 
errors f z = 1, f z = 2, f x = 1 and f x = 2, 

^El > 1 - Vz{l) - T} Z {2) - 77x(l) - r)x(2) = 0.792. 

(11) 



It is therefore possible to improve the estimate of en- 
tanglement capability originally reported in [2] with- 
out adding any new data. 

It may also be interesting to apply the same anal- 
ysis to the reverse operation, that converts the Bell 
states into local XZ eigenstates. The fidelity F xz of 
this disentangling operation is given by the process 
matrix element sum 

F xz = F qp + X(2,0)(2,0) +X(0,1)(0,1) + X(2,l)(2,l) ! (I 2 ) 

and the minimal fidelity defined by the measurement 
results is 

F XZ >1- rizQ) ~ Vz(3) - Vx(2) - rjx{3) = 0.806. 

(13) 

The minimal fidelity of disentangling the Bell states is 
thus higher than the minimal fidelity of entanglement 
generation by the time-reversed process. 

Finally, it is also possible to consider the errors 
in entanglement generation from YY inputs. The 
outputs are then maximally entangled states with the 
stabilizers //, ZY, YX, and XZ. Therefore, the 
fidelity F g2 is given by 

^E2 = ^gp + X(l,3)(l,3) +X(3,2)(3,2) +X(2,1)(2,1)- (14) 

In this case, all of the stabilizers represent errors that 
show up in both the Z and the X measurements. 
Therefore, the fidelity F xz can be as low as the min- 
imal process fidelity, 

F E2 > Fz + F x - 1 = 0.720. (15) 

Thus the measurement data provided by the comple- 
mentary Z and X operations provides only a very 
rough estimate of the operation on Y inputs. 

Of course, the worst case noise estimate given in ta- 
ble Q] is a rather extreme and unlikely interpretation 
of the measurement data. It may therefore be inter- 
esting to compare it with a more realistic estimate 
based on the assumption that the errors in Z and 
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Table 4: "Worst case" estimate of process matrix elements for the experimental data from 2 shown in table 



XtfzJxUfzJx) 


fx = 


fx = l 


fx = 2 


/* = 3 


Sum 


fz = 


0.720 


0.034 


0.071 


0.028 


0.853 


f z = l 


0.051 











0.051 


f z = 2 


0.052 











0.052 


/, = 3 


0.044 











0.044 


Sum 


0.867 


0.034 


0.071 


0.028 


1 



X are uncorrelated @]. For this purpose, it is useful 
to consider the definition of the average fidelity of a 
quantum process, 

F av . = /d$($ | lllE(\ \)U Q | $), (16) 



where the integral over $ represents a uniform av- 
erage over all possible states of the (/-dimensional 
Hilbert space. It has been shown |S] that this av- 
erage fidelity is related to the process fidelity by 

_ F qp d+l 

Fav - - "^Ti - - (17) 

This relation can be explained quite intuitively by 
using the error analysis above. As eq.(@J indicates, 
the fidelities for one input basis of d orthogonal states 
are given by a sum of the process fidelity F qp and the 
diagonal process matrix elements of d—1 errors, out of 
a total number of d 2 — l = (d+l)(d—l) possible errors. 
This means that the probability that the operation on 
any given input state | $) is insensitive to any given 
error is (d - l)/(d 2 - 1) = l/(d + 1). The average 
over all possible input states | $) makes this relation 
exact: the contribution of l/(d+ 1) to F a v. in eq.JTJJ 
simply represents the probability that the errors in 
the quantum operation result only in an unobservable 
phase change. 

This interpretation of eq. I|17|) indicates that we can 
consider the complementary fidelities Fx and Fz as 
representative fidelities contributing to the average 
fidelity Fav.- If the complementary fidelities are not 
too different from each other, the most realistic as- 
sumption seems to be that the average fidelity -Fav. is 
close to the average of the complementary fidelities. 
An estimate of the most likely value of the process 
fidelity F qp can then be obtained from 



Fqp (est. 



1 

1+ d 



(Fz_ 



Fx 



(18) 



Since this estimate implicitly assumes a rather uni- 
form distribution of errors, it is natural to apply 



the same assumption to obtain an estimate for the 
remaining diagonal elements of the process matrix. 
This means that the joint probabilities of errors given 
by X(fx,fx)(f*,fx) should be proportional to the prod- 
uct of the error probabilities r)z(f z ) and rjx(fx) for 
fz, fx 7^ 0. With this assumption of uncorrelated er- 
rors, the remaining process matrix elements read 



X(/„/.)c/.,/.)(ffltO 



d- 1 / 1 




i 


d \1-Fz~* 


1 


-F x 


X(/,,o)(/„,o)(est.) ~ 






fd+1 d-l 


n 


~F X 


V 2d 2d 




-F z 


X(o,/«)(o,/.)(est.) » 






fd+l d-l 


n 


-F z 



Vz(fz)vx(fx) 
Vz(fz) 

Vx(f x ). (19) 



V 2d 2d V 1 - F x 

This estimate may be the most realistic description 
for quantum processes with very similar complemen- 
tary fidelities Fz and Fx ■ If the values of Fz and Fx 
are very different, it is possible that the estimates for 
X(o,/ x )(o,/ x ) or for X(/„o)(/„o) become negative, espe- 
cially if the higher fidelity is close to one. Although 
such errors can be corrected by restricting all diag- 
onal elements to positive values, it would be more 
natural to use a different noise models for such cases. 

In the case of our controlled-NOT gate, the similar 
values of Fz and Fx indicate that the statistical noise 
model may be appropriate. The corresponding dis- 
tribution of process matrix elements is shown in table 
\5\ All diagonal elements representing errors now have 
very similar values, ranging from a minimum of 0.59 
% for (0, 3) to a maximum of 1.98 % for (2, 2). The es- 
timated process fidelity of F qp — 0.825 is only a little 
bit lower than the maximal possible process fidelity 
of Min{Fz, Fx} = 0.853, and much higher than the 
lower limit of Fz + Fz — 1 = 0.720. Statistical consid- 
erations thus indicate that the actual process fidelity 
of the device is likely to significantly exceed the min- 
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Table 5: Statistical estimate of the process matrix elements for the experimental data from 2. shown in table 
[3 The assumption of uncorrelated errors results in a uniform distribution of errors over all possibilities. 





fx = 


fx = l 


fx = 2 


/* = 3 


Sum 


fz = 


0.825 


0.0072 


0.0150 


0.0059 


0.853 


f z = l 


0.0146 


0.0093 


0.0194 


0.0077 


0.051 


fz = 2 


0.0149 


0.0095 


0.0198 


0.0078 


0.052 


A = 3 


0.0126 


0.0080 


0.0168 


0.0066 


0.044 


Sum 


0.867 


0.034 


0.071 


0.028 


1 



imum assumed in the "worst case" estimate shown in 
table H 

We can also derive estimates for other operations 
from the statistical noise model, corresponding to 
their most likely values. The results derived from 
eqs. (0 HDnHHl read 

F zx « 0.874 F m « 0.850 

F X z ~ 0.857 F E2 « 0.859. (20) 

Naturally, all of these fidelities are now close to the 
average fidelity of F av . = (F z + F x )/2 = 0.86. Inter- 
estingly, the estimate for the fidelity F^2 of entangle- 
ment generation from YY inputs is now higher than 
the estimate for the fidelity of entanglement gen- 
eration from XZ inputs. This change illustrates the 
fundamental difference between minimal values and 
likely values in the fidelity estimates. In fact, the two 
results show that the available data allows a far more 
precise estimate of F^ than of -Fgo- 

5 Conclusions 

The errors observed in A-qubit operations can be 
characterized in terms of the bit flip errors observed 
in the complementary operations resulting in Z and 
in X output states. It is therefore useful to char- 
acterize a device by first measuring the 2d = 2 N+1 
output fidelities of these two complementary opera- 
tions. Since the only process able to perform both 
operations with a fidelity of 1 is the intended process 
Uq, it is possible to derive upper and lower bounds 
for the process fidelity from these two measurements. 
For high fidelity processes, this estimate is sufficient 
to confirm the successful implementation of a multi 
qubit gate. 

For noisy processes such as the quantum 
controlled-NOT analyzed above, the details results 
for the fidelities and errors allow estimates of the pro- 
cess matrix elements corresponding to various noise 



models. Even though these error models are not as 
precise as the estimates obtained from full quantum 
tomography, it is remarkable that such a detailed 
analysis is possible using only a small fraction of 
the d 4 = 16 N measurement probabilities required for 
complete quantum process tomography. The evalua- 
tion of complementary operations is therefore a par- 
ticularly efficient method for the characterization of 
multi-qubit quantum devices. 
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